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AERONAUTIC  SYMBOLS 


1.  FUNDAMENTAL  AND  DERIVED  UNITS 


Symbol 

Metric 

■ 

English 

Unit 

Abbrevia¬ 

tion 

Unit 

Abbrevia¬ 

tion 

Lenfftli _ 

1 

meter _ _ _ 

m 

foot  (or  mUe) _ _ _ _ 

ft  (or mi) 

Time _ _  _ 

t 

second _ 

s 

second  (or  hour).  _  __ 

sec  (or  hr) 

Force _ 

F 

weight  of  1  kilogram _ ; 

kg 

weight  of  1  pound _ 

lb 

Power _ ; _ _ 

P 

horsepower  (metric) ____ 

horsepower _ 

hp 

n  1 

V 

/kilometers  per  hour _ 

kph 

miles  per  hour _ _ 

mph 

[meters  per  second _ 

mps 

feet  per  second _ _ 

fps 

2.  GENERAL  SYMBOLS 

Weight==jn^  ^  ^  _  _  v  Kinematic  viscosity 

Standard  acceleration  of  gravity=9.80665-m/a*  p  Density  (mass  per  unit  volume) 
or  32.1740  ft/sec*  Standard  density  of  dry  air,  0.12497  kg-m“*-s®  at  15°  C 

,,  _W'  -  and. 760  mm;  or  0.002378  Ib-ft'^  sec^® 

mass— ^  Specific,  weight  of  “standard”  air,  1.2265  kg/m®  or 

Moment  of  me:rtia=mP.  (Indicate  ^s  of  0.07661  Ib/cu  ft  _ 

radius  of  gyration  k  by  proper  subscript.) 

'Coefficient  of  viscosity 

3.  AERODYNAMIC  SYMBOLS 


Area 

Area  of  wing 
Gap 
\Span 
Chord 

Aspect  ratio,  ' 

True  air  speed 
Dynamic  pressure,  F® 

Lift,  absolute  coefficient 

Drag,  absolute  coefficient 

Profile  drag,  absolute  coefficient 

Induced  drag,  absolute  coefficient  Cb.=^ 

‘go 

Parasite  drag,  absolute  coefficient  Cdp=—^ 


Angle  of  setting  of  wings  (relative  to  thrust  line) 
Angle  of  stabiliaer  setting  (relative  to  thrust 
line)  -  : 

Resultant  moment 
Resultant  angular  velocity 

Reynolds  nmnber,  p~  where  Z  is  a  linear  dimen- 

sion  (e.g.,  for  an  airfoil  of  1 .0  ft  chord,  100  mph, 
standard  pressure  at  15°  C,  the  corresponding 
Reynolds  nmhber  is  935,400;  or  for  an  airfoil 
of  1.0  m  chord,  lOO  mps,  the  corresponding 
Reynolds  munber  is  6,865,000) 

Angle  of  attack 
Angle  of  downwa^ 

Angle  of  attack,  infinite  aspect  ratio 
Angle  of  attack,  induced 
Angle  of  attack,  absolute  (measmed  from  zero- 
lift  position)  , 

Fhght-path  angle 


Cross-wind  force,  absolute  coefficient  Gc=- 


tech  ubrary  kafb,  nm 

■llllllll 


014315*! 


REPORT  No.  882 


FREQUENCY-RESPONSE  METHOD  FOR  DETERMINATION 
OF  DYNAMIC  STABILITY  CHARACTERISTICS  OF 
AIRPLANES  WITH  AUTOMATIC  CONTROLS 

By  HARRY  GREENBERG 


Langley  Memorial  Aeronautical  Laboratory 
Langley  Field,  Va. 


National  Advisory  Committee  for  Aeronautics 

Headquarters,  1724  F  Street  NW,  Washington  25,  D.  C. 

Created  by  act  of  Congress  approved  March  3,  1915,  for  the  supervision  and  direction  of  the  scientific  study 
of  the  problems  of  flight  (U.  S.  Code,  title  49,  sec.  241).  Its  membership  was  increased  to  15  by  act  appi'oved 
March  2,  1929.  The  members  are  appointed  by  the  President,  and  serve  as  such  without  compensation. 

Jerome  C.  Hunsaker,  Sc.  D.,  Cambridge,  Mass.,  Chairman 

Alexander  Wetmore,  Sc.  D.,  Secretary,  Smithsonian  Institution,  Vice  Chairman 


Hon.  John  R.  Alison,  Assistant  Secretary  of  Commerce. 

Vannevar  Bush,  Sc.  D.,  Chairman,  Research  and  Development 
Board,  Department  of  Rational  Defense. 

Edward  U.  Condon,  Ph.  D.,  Director,  National  Bureau  of 
Standards. 

Donald  B.  Duncan,  Vice  Admiral,  Deputy  Chief  of  Naval 
Operations  (Air). 

R.  M.  Hazen,  B.  S.,  Chief  Engineer,  Allison  Division,  General 
Motors  Corp. 

William  Littlewood,  M.  E.,  Vice  President,  Engineering, 
American  Airlines  System. 

Theodore  C.  Lonnquest,  Rear  Admiral,  Assistant  Chief  for 
Research  and  Development,  Bureau  of  Aeronautics,  Navy 
Department. 


Edward  M.  Powers,  Major  General,  United  States  Air  Eorce, 
Deputy  Chief  of  Staff,  Materiel. 

Arthur  E.  Raymond,  M.  S.,  Vice  President,  Engineering, 
Douglas  Aircraft  Co. 

Francis  W.  Reichelderfer,  Sc.  D.,  Chief,  United  States 
Weather  Bureau. 

Carl  Si’aatz,  General,  Chief  of  Staff,  United  States  Air  Force. 
Orville  Wright,  Sc.  D.,  Dayton,  Ohio. 

Theodore  P.  Wright,  Sc.  D.,  Administrator  of  Civil  Aero¬ 
nautics,  Department  of  Commerce. 


Hugh  L.  Dryden,  Ph.  D.,  Director  of  Aeronautical  Research 
John  W.  Crowley,  Jr.,  B.  S.,  Associate  Director  of  Aeronautical  Research 


John  F.  Victory,  LL.M.,  Executive  Secretary 
E.  H.  Chamberlin,  Executive  Officer 


Henry  J.  E.  Reid,  Sc.  D.,  Director,  Langley  Memorial  Aeronautical  Laboratory,  Langley  Field,  Va. 

Smith  J.  DeFrance,  B.  S.,  Director  Ames  Aeronautical  Laboratory,  Moffett  Field,  Calif. 

Edward  R.  Sharp,  LL.  B.,  Director,  Flight  Propulsion  Research  Laboratory,  Cleveland  Airport,  Cleveland,  Ohio 


TECHNICAL  COMMITTEES 


.\Ell01>YNAMiC.S 

Power  Plants  for  Aircraft 

.  V I KCR  A  FT  Co  .N  STR  UC’ITO  N 


Operating  Probi.ems 
Self-Propelled  Guided  Mlssiles 
Industry  Con.sulting 


Coordination  of  Research  Needs  of  Military  and  Civil  Aviation 
Preparation  of  Research  Programs 
Allocation  of  Problems 
Prevention  of  Duplication 
Consideration  of  Inventions 


Langley  Memorial  .\euonauttcal  Laboratory,  Ames  .Aeronautical  Laboratory, 

Langley  Field,  Va.  Moffett  Field,  Calif. 

Flight  Propulsion  Research  Laboratory, 

Cleveland  Airport,  Cleveland,  Ohio 

Conduct,  under  unified  control,  for  all  agencies,  of  scienlifie  research  on  the  fundamental  problems  of  fight 


ll 


Office  of  Aeronautical  Intelligence, 

Washington,  D.  C. 

Collection,  classification,  compilation,  and  disset lination  of  scientific  and  technical  information  on  aeronautics 


REPORT  No.  882 


FREQUENCY-RESPONSE  METHOD  FOR  DETERMINATION  OF  DYNAMIC  STABILITY 
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By  Harry  Greenberg 


SUMMARY 

A  frequency-response  method  for  determining  the  critical 
control-gearing  and  hunting  oscillations  of  airplanes  with 
automatic  pilots  is  presented.  The  method  is  graphical  and 
has  several  advantages  over  the  standard  numerical  procedure 
based  on  Routh’s  discriminant.  The  chief  advantage  of  the 
method  is  that  direct  use  can  he  made  of  the  measured  response 
characteristics  of  the  automatic  pilot.  This  feature  is  especially 
useful  in  determining  the  existence,  amplitude,  and  frequency 
of  the  hunting  oscillations  that  may  be  present  when  the  auto¬ 
matic  pilot  has  nonlinear  dynamic  characteristics . 

Several  examples  are  worked  out  to  illustrate  the  application 
of  the  frequency-response  method  in  determining  the  effect  of 
automatic-pilot  lag  or  lead  on  critical  control  gearing  and  in 
determining  the  amplitude  and  frequency  of  hunting.  It  is 
shown,  that  the  method  may  be  applied  to  the  ca.se  of  a  control 
geared  to  airplane  motions  about  two  axes. 

INTRODUCTION 

The  increased  use  of  automatic  control  on  aircraft  (espe¬ 
cially  on  pilotless  aircraft)  lias  focused  attention  on  the 
proper  design  of  the  aircraft  and  associated  control  systems 
with  a  view  toward  obtaining  satisfactory  dynamic  stability. 
The  control  system  (autopilot)  consists  of  the  gyro,  phase- 
shifting  device,  and  servomotor.  Factoi's  in  the  control 
system  that  determine  stability  characteristics  are  control 
gearing,  lag  in  the  servomotor,  and  lead  in  the  phase- 
shifting  device.  The  stability  characteristies  of  the  air¬ 
plane  depend  on  airplane  configuration  and  mass  distribution. 

The  purpose  of  the  present  report  is  to  give  a  method  for 
analyzing  the  dynamic  stability  of  an  airplane  with  auto¬ 
matic  control  which  makes  direct  use  of  the  observed  dynamic 
characteristics  of  the  automatic  pilot  (herein  called  auto¬ 
pilot).  The  method  separates  the  characteristics  of  the 
autopilot  from  those  of  the  airplane.  It,  therefore,  easily 
reveals  the  effects  of  modifications  to  the  autopilot  such  as 
adding  lead.  The  procedure  is  largely  graphical.  The 
method  is  similar  in  certain  respects  to  that  of  Nyquist 
(reference  1),  which  was  devised  for  electronic  circuits  but 
which  has  also  been  applied  by  some  workers  to  the  design 
of  servomotors. 

The  frequency-response  method  was  previously  applied  by 


Jones  (reference  2)  to  calculate  the  hunting  produced  in 
airplanes  with  “flicker,”  that  is,  on-off  control.  If  the  servo¬ 
motor  (herein  called  servo)  has  a  linear  lag  characteristic 
(lag  independent  of  amplitude),  then  there  will  usually  be  an 
upper  limit  to  the  control  gearing  above  which  the  airplane 
will  be  unstable.  The  frequency-response  method  deter¬ 
mines  this  critical  control  gearing  and  the  corresponding 
frequency  for  any  type  of  lag  or  lead  in  the  autopilot.  It 
also  indicates  the  changes  to  be  made  in  the  servo  which  will 
improve  the  stability  of  the  airplane. 

The  utility  of  the  frequency-response  method  is  most  ap¬ 
parent  in  the  ease  of  an  actual  servo  with  nonlinear  lag 
characteristies.  This  mc'thod  allows  the  use  of  measured 
characteristics  of  the  servo  that  might  be  inconvenient  to 
represent  by  a  mathematical  formula. 

TERMINOLOGY  AND  SYMBOLS 

The  word  “lead”  is  used  in  this  report  in  two  ways.  It  is 
the  phase  angle  of  control-surface  dcfleetion  6  ahead  of  air¬ 
plane  deflection,  say  angle  of  pitch  9,  and  it  is  also  used  to 
indicate  a  device  that  causes  an  increase  in  the.  phase  angle 
of  lead  of  5  ahead  of  9.  Those  phase-shifting  devices  are 
herein  called  first-derivative  lead  and  second-derivative  load. 

The  expression  “linear  autopilot”  is  used  to  indicate  a  servo 
that  is  acted  upon  by  a  force  proportional  to  the  input  signal 
(airplane  deflection)  and  resisted  by  a  force  proportional  to 
the  displacement  and  velocity  of  the  output  (control-surface 
motion). 

The  function  of  the  automatic  pilot  is  to  apply  a  eorreetive 
control  deflection  in  response  to  any  deflection  of  the  air¬ 
plane.  The  “control  gearing”  is  the  ratio  of  the  applied  con¬ 
trol  deflection  to  the  airplane  deflection  for  very  slow 
deflections  (static  condition).  There  may  exist  a  “critical 
(iontrol  gearing,”  which  in  nearly  all  practical  eases  is  an 
upper  limit  beyond  wliich  dynamic  instability  occurs. 

The  following  symbols  are  used: 
a  ratio  of  control  deflection  due  to  second  derivative  of 
airplane  displacement  to  that  due  to  airplane 
displacement 

r  ratio  of  control  deflection  due  to  first  derivative  of 
au’plane  displacement  to  that  due  to  airplane 
displacement 
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R  amplitude  of  airplane  oscillation  produced  by  unit 
amplitude  oscillation  of  control  surface 
Sr  amplitude  of  control-surface  oscillation  required  to 
maintain  unit  amplitude  of  airplane  oscillations  (1/fi) 
Sp  amplitude  of  control-surface  oscillation  produced  by 
autopilot  in  response  to  airplane  oscillation,  divided 
by  amplitude  of  airplane  oscillation 
Spj  5p  for  unit  control  gearing 
kcr  critical  value  of  control  gearing 
ki  spring  stiffness  factor  of  servo  and  control  surface 

ki  factor  relating  airplane  deflection  and  force  on  servo 

k  control  gearing,  ratio  of  S  to  6  in  static  condition 

(k.lkO 

c  viscous  damping  factor  of  servo  and  control  surface 
I  viscous  lag  factor  which  depends  on  c  and  ki  (c/k,) 

m  mass  of  movable  part  of  servo  and  control  surface 

t  time,  seconds 

T  period  of  oscillation,  seconds 

6  angle  of  control-surface  deflection  in  dii'ection  f.o 

reduce  angle  of  pitch  of  airplane 
tp  phase  angle  of  lead  of  &  ahead  of  6  when  oscillating 
au'plane  forces  control  surface  to  oscillate 
ir  phase  angle  of  lead  of  S  ahead  of  6  when  oscillating 
control  surface  forces  airplane  to  oscillate 
8  angle  of  pitch  of  airplane 

0  signal  fed  into  servo 

<!>  angle  of  bank  of  airplane 

i/'  angle  of  yaw  of  airplane 

$  gyro  displacement  produced  by  airplane  deflection 

w  angular  frequency,  radians  per  second  (2x  XFrequency) 

oi„  angular  frequency  of  servo  and  airplane  when  k=kc.r, 
radians  per  second 

natural  angular  frequency  of  servo,  radians  per  second 
r  angle  of  tilt  of  gyro  axis  from  X-axis  of  airplan(^ 

D  differential  operator  {djdt) 

yubscript; 

max  maximum 

BASIC  PRINCIPLES  OF  FREQUENCY-RESPONSE  METHOD 

The  response  of  the  aiiplane  to  a  sinusoidal  control  motion 
may  be  measured  in  flight  or  may  be  computed  from  the 
equations  of  motion  of  the  airplane.  (See  appendix  A.) 
The  airplane  response  is  sinusoidal  and  of  the  same  frequency 
as  the  control  motion  and  has  an  amplitude  and  phase  that 
depends  on  this  frequency.  If  the  control  motion  is  given  by 

5=sin  co< 

then  the  response  in  pitch,  for  example,  may  be  expressed  as 
d=R  sin  (coi—  e,) 

where  R  and  er  both  depend  on  co,  and  R  is  the  ratio  of  the 
amplitude  of  6  to  the  amplitude  of  6  when  the  control  surface 
is  forcing  the  airplane  to  oscillate.  Changing  the  amplitude 
of  S  produces  a  proportionate  change  in  the  amplitude  of  9 
but  does  not  affect  the  phase  angle  e^.  The  reciprocal  of  R 


gives  the  amplitude  of  5  required  to  sustain  unit  amplitude 
oscillation  in  8  and  is  denoted  by  Sr]  that  is,  if 

S=Sr  sin  wt 

then  the  response  will  be 

0=sin  (wt — e,) 

A  plot  of  Sr  and  against  co  is  then  made  and  is  referred  to 
as  the  frequency  response  of  the  airplane. 

The  frequency  response  of  the  autopilot  is  obtained  by 
measuring  or  computing  the  response  of  the  control  surface 
to  a  steady  oscillation  of  the  airplane.  The  response  may 
be  calculated  if  the  dynamic  constants  of  the  autopilot  are 
known,  but  it  is  usually  easier  to  measure  the  response  by 
oscillating  the  autopilot  in  the  laboratory.  These  calcula¬ 
tions  or  measurements  give  the  response  to 

0=sin  wt 

and  this  response  is  expressed  by 

S  =  Sp^  sin  (cof-h  ep) 

where  Sp^  is  the  amplitude  of  S  for  unit  control  gearing  and 
ep  is  the  phase  angle  of  lead  of  S  ahead  of  9  when  the  airplane 
is  forcing  the  control  surface  to  oscillate. 

If  at  some  value  of  the  angular  frequency  w  the  condition 

€p  ^  €r 

exists,  then  the  airplane  connected  to  the  autopilot  will 
oscillate  continuously  at  this  frequency  at  constant  ampli¬ 
tude,  provided  the  control  gearing  is  made  equal  to  Sr/Sp^  at 
this  value  of  a>;  that  is,  the  critical  control  gearing  is  given  by 


at  ft)  where  er=tp.  This  condition  of  neutral  stability  is 
shown  schematically  in  figure  1.  An  example  of  a  solution 
for  kcr  is  shown  in  figure  2.  In  this  case  Wct=&.2  and 

5  0  33 

kcr=}r  at  this  value  of  co;  therefore  A:c,=fr^  =  0.36. 

u .  y  I 

Assumec/ 


Figure  1. — Response  of  airplane  to  sinusoidal  control  motion  and  response  of  control  to 
sinusoidal  airplane  motion.  If  €p  =  <r  find  k=kcT—^}  the  autopilot  will  sustain  the 

dpj 


oscillations. 
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Figure  2.— Illustration  of  froquoncy-response  method  of  determining  critical  control  gearing, 
au,=6,2;  kcr=0.M.  (Dotted  line  indicates  whore  <p=«,.) 


In  almost  all  cases  control  gearings  greater  than  kc,  result 
in  instability  and  control  gearings  less  than  kc,  provide 
stability.  Thus,  k,,  is  usually  an  upper  limit  to  the  permis¬ 
sible  control  gearing.  More  generally,  the  stable  side  of 
the  boundary  is  determined  by  the  relative  slopes  of  the  t. 


and  curves  at  their  intersection. 


If  ^  is  greater  than 


the  system  is  stable  below  kc,  and  vice  versa.  A  proof 

of  this  rule  is  given  in  appendix  B. 

In  the  preparation  of  figure  2  and  the  examples  that  follow 
(figs.  3  to  8)  a  particular  airplane  configuration  was  assumed 
for  the  computation  of  8,  and  t,.  (See  appendix  A  for 
methods.)  This  hypothetical  pilotless  airplane  has  four 
fins  90°  apart,  a  horizontal  pair  for  pull-ups  and  a  vertical 
pair  for  turns.  The  fins  constitute  the  only  lifting  surfaces 
of  the  airplane  and  are  equipped  with  ailerons  connected  to 
a  roll  gyro  and  servo  that  maintain  the  airplane  in  a  position 
in  which  one  pair  of  fins  is  always  horizontal.  The  pitch 
and  yaw  control  surfaces  are  trailing-edge  flaps  on  the  fins 
and  are  connected  in  pairs,  one  pair  on  the  horizontal  fins 
connected  to  a  pitch  gyro  and  servo  and  the  other  pair  on 
the  vertical  fins  connected  to  a  yaw  gyro  and  servo.  It  is 
assumed  that  the  pitch  and  yaw  motions  do  not  interact 
with  each  other  or  with  any  accidental  rolling  motion  that 
826251 — 49 - 2 


Figure  3.— ElTeet  of  lag  in  the  servo  on  tlie  frequenoy  response  of  a  critically  damijod  linear 

2 

servo.  N  alural  frequency  of  servo 

may  exist.  Figures  2  to  6  show  the  elevator-pitch  stability. 
Rudder-yaw  stability  may  be  expected  to  be  very  nearly 
the  same  as  elevator-pitch  stability  because  of  the  configura¬ 
tion  of  the  fins. 

LINEAR  AUTOPILOT 

The  motion  of  a  control  surface  actuated  by  a  simple  linear 
autopilot  is  similar  to  that  of  a  mass-spring-dashpot  and  may 
be  represented  by  the  equation 

mD'^8-\-cDS-[-ki8=k26  (1) 

where  each  of  the  factors  m,  c,  and  ki  is  the  sum  of  two  parts, 
one  part  due  to  the  servo  and  one  part  due  to  the  control 
surface.  For  example,  ki  is  the  sum  of  the  spring  constant 
of  the  servo  and  the  aerodynamic  hinge-moment  constant 
of  the  control  surface.  The  term  k26  represents  a  force 
proportional  to  the  deviation  in  angle  of  pitch  6  applied 
without  lag. 

Actually,  the  effect  of  control-cable  flexibility  is  to  intro¬ 
duce  an  additional  degree  of  freedom  into  the  system.  This 
effect  is  probably  small  and  may  be  taken  into  account 
either  analytically  or  by  flight  measurements  as  mentioned 
in  the  section  “Nonlinear  Autopilot.’' 
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Equation  (1)  may  be  expressed  as 

where  k^jki  is  the  control  gearing  k.  The  natural  frequency 
of  the  servo  is  given  by 


and  the  damping  depends  on  the  quantity  With  these 

substitutions,  equation  (2)  becomes 

5(\-\-lD+~^=^kB 

\  “n  / 

The  term  ID  represents  viscous  damping  and  the  term 
D'^lwn  represents  inertia  reaction. 

When  the  autopilot  is  oscillated  in  pitch,  the  control 
response  S  lags  behind  6.  The  lag  usually  decreases  as  the 
lag  factor  I  decreases  and  as  the  natural  frequency  co„ 


Figure  5.— Effect  of  lag  and  lead  on  the  critical  control  gearing  with  a  critically  damped 


linear  servo.  Airplane  characteristics  of  figure  2. 


l+rD 


increases.  It  is  convenient  to  choose  equal  to  2/Z,  a  value 
that  gives  critical  damping.  With  this  choice,  the  phase 
angle  and  the  relative  amplitude  for  unit  control  gearing 
of  the  response  are  presented  in  figure  3  for  four  values 
of  1.  Values  of  5r  and  e,  for  the  assumed  airplane  are  also 
given  so  that  kcr  can  be  obtained  for  each  value  of  1. 

If  first-derivative  lead  is  added  to  the  system,  the  equation 
of  motion  relating  5  and  6  becomes 

5  (l  +  lD+~^=k{l+rD)d 
\  / 


Figure  4  shows  and  tp  for  the  same  values  of  I  and  as 

figure  3,  but  with  Comparison  shows  that  the  first- 

derivative  lead  produces  an  angle  of  lead  that  approaches 
90°  at  high  values  of  the  angular  frequency  co.  This  amount 
of  lead  is  enough  to  override  the  lag  for  moderate  values  of  I 
at  low  frequencies  but  not  at  high  frequencies. 

Values  of  kcr  obtained  from  figures  3  and  4  are  shown  in 
figure  5  plotted  against  the  viscous  lag  factor  1.  Considerable 
increase  in  the  stable  range  of  control  gearing  is  evident  for 
small  values  of  I  with  the  first-derivative  lead.  Note  how 
this  improvement  in  stability  diminishes  at  large  values  of  1. 

If  both  first-derivative  lead  and  second-derivative  lead 
are.  added,  the  response  of  the  control  is  given  by  the  equation 


5=k 


2)2“  ® 

1-f  ZD-f — 2 

<^n 


NONLINEAR  AUTOPILOT 

Even  if  the  servo  is  constructed  to  give  a  response  pro¬ 
portional  to  a  steady  disturbance,  its  response  to  an  oscillating 
disturbance  is  not  proportional  in  practice  to  the  amplitude 
of  the  disturbance  because  of  the  nonlinear  dynamic  character¬ 
istics  of  the  servo.  The  lag  is  then  a  function  of  both  am¬ 
plitude  and  frequency.  By  finding  experimentally  the 
frequency  response  of  the  autopilot  for  a  number  of  amplitudes 
and  for  a  given  control  gearing,  it  is  possible  to  determine  the 
amplitude  and  frequency  at  which  the  airplane  will  oscillate 
when  coupled  to  the  autopilot  at  that  particular  value  of 
control  gearing.  The  condition  for  steady  oscillations 
(hunting)  is  that  dp=8r  and  ep=  at  some  frequency.  The 
frequency  and  amplitude  of  the  steady  oscillations  of  the 
airplane  are  the  values  of  co  and  dmax  at  which  3^=5^  and 


case  of  the  hunting  that  may  exist  with  a  nonlinear  servo 
the  oscillations  are  said  to  be  stable  if,  after  a  disturbance 
from  their  steady  value,  they  tend  to  return  to  that  steady 
value.  The  criterion  for  the  stability  of  the  steady  oscilla¬ 
tions  follows  from  that  for  the  linear  autopilot  and  is 

d  — 

djer—ep) _ ^ 

dco  dOf/mx 

The  solution  obtained  is  not  exact  because  the  autopilot 
response  to  sinusoidal  airplane  motion  is  not  a  pure  sine  wave 
but  is  distorted;  however,  the  existence  and  approximate 
amplitude  and  frequency  of  hunting  oscillations  of  the  funda¬ 
mental  can  be  determined  by  the  frequency-response  method. 
The  values  of  8p  and  used  when  the  response  is  not  sinusoi¬ 
dal  are  the  values  of  the  equivalent  sinusoidal  reponse  which 


ep=  ij. 
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Figure  8.— Plot  of  5p/Sr  and  w«r  against  dmax  obtained  from  figure  7,  showing  the  amplitude 
and  frequency  of  the  hunting. 

supplies  the  same  energy  per  cycle  and  the  same  impulse 
per  half  cycle  as  the  observed  response.  This  condition  leads 
to  the  following  expressions  for  and  5^: 

,  -1^ 

«P=tan  2 

where 


The  value  of  5  in  these  integrals  is  the  measured  control- 
surface  deflection  in  response  to  an  oscillating  motion  of  the 
airplane  represented  by  6=dmax  sin  ut.  It  turns  out  that  B 


is  exactly  equal  to  the  cosine  coefficient  of  the  fundamental 
frequency  component  of  the  Fourier  series  approximation 
to  S;  however,  A  is  not  in  general  exactly  equal  to  the  first 
sine  coefficient  of  the  Fourier  series.  Ordinarily  the  oscilla¬ 
tion  in  pitch  is  performed  in  the  laboratory  by  mounting  the 
autopilot  on  a  table  that  can  be  made  to  oscillate. 

A  more  exact  method  of  determining  the  autopilot  charac¬ 
teristics,  which  includes  the  effects  of  control-system  flex¬ 
ibility  and  hinge  moment,  is  to  measure  the  response  of  the 
control  snrface  and  airplane  in  flight  to  a  sinusoidal  input 
signal  to  the  servo.  In  order  to  obtain  the  frequency  response 
of  the  control  system  to  sinnsoidal  input  signal,  it  is  necessary 
to  divide  the  vector  ontpnt  S  by  the  vector  sum  of  input 
signal  dne  to  6  and  inpnt  signal  applied  directly.  Thus,  the 
dynamic  characteristics  of  the  airplane  and  control  system 
can  be  obtained  at  the  same  time.  If  the  aircraft  is  not 
available  for  flight  tests,  one  can  “synthesize”  the  complete 
servo-control-system  response  from  gronnd  measurements  of 
the  servo  and  estimates  of  control-system  flexibility,  inertia, 
and  aerodynamic  hinge  moments. 

A  typical  application  of  the  method  for  determining  ampli¬ 
tude  and  freqnency  of  hnnting  of  an  aircraft  and  antopilot 
with  stops  in  the  control  system  is  shown  in  fignres  7  and  8. 
The  airplane  nsed  in  this  example  is  slightly  different  from 
that  considered  previously  and,  therefore,  the  Sr  and  e,  curves 
are  slightly  different  from  those  in  figures  2  to  6.  Figure  7 
contains  the  phase  and  amplitude  of  the  experimental 
response  of  the  autopilot  plotted  against  frequency  for  a 
series  of  amplitudes  of  the  input  and  also  the  corresponding 
calculated  curves  for  the  aircraft.  For  each  amplitude  the 
values  of  dp  and  Sr  are  obtained  at  the  frequency  where 
€p=fr-  The  ratio  Sp/Sr  and  the  value  of  (where  ep=fr)  are 
plotted  against  amplitude  of  pitch  Bmax  in  figure  8.  The  ampli- 

tude  and  frequency  of  hunting  occur  where  ^=1-  In  this 

0  r 

example  the  amplitude  of  the  hunting  is  about  2.5°  and  the 
frequency  is  18.7  radians  per  second. 

Nonlinearity  in  the  autopilot  usually  causes  an  increase  of 
lag  with  amplitude.  This  increase  may  result  in  a  reversal  of 
the  favorable  effect  of  first-derivative  lead  as  the  lead  factor 
r  is  indefinitely  increased  and  suggests  an  optimum  value  for 
r,  a  conclusion  that  is  supported  by  flight  experience.  This 
result  may  be  traced  to  the  counteracting  effects  of  lead  and 
amplitude  introduced  by  the  first-derivative  lead.  The  first- 
derivative  lead  shifts  the  phase  of  the  signal  fed  into  the  servo 
so  that  the  signal  leads  B  by  an  amount  tan”i  rco  and  increases 
the  amplitude  of  the  signal  by  the  factor  The 

change  of  lead  with  r  may  be  stated  symbolically  as  follows; 
If  the  dependence  of  on  B  and  co  without  first-derivative 
lead  is  expressed  by 

ip=JiB,  w) 

then  the  addition  of  first-derivative  lead  changes  the  lead  to 
ep=/(0,  a))-)-tan“Vo  (4) 
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In  order  to  apply  the  frequency-response  method  in  this  ease 
it  is  necessary  only  to  calculate  the  response  in  bank  and 
yaw  of  the  airplane  to  a  sinusoidal  control  motion  and  to 
combine  them  according  to  the  preceding  formula  to  obtain 
This  formula  yields  the  5,  and  curves.  The  other  pair 

of  curves  required — the  Sp  and  curves — are  obtained  in 

exactly  the  same  way  as  before,  namely,  by  oscillating  the 
gyro  about  its  sensitive  axis.  The  critical  control  gearing  so 
determined  will  be  the  critieal  ratio  between  5  and  '!>. 


where 

Differentiating  equation  (4)  with  respect  to  r  gives 

dip  CO 

^dd  vr+^^"'"r+7i? 

If  r  IS  very  small,  the  first  term  on  the  right-hand  side  of 
equation  (5)  is  negligible,  so  that  ^  is  positive.  The  addi¬ 
tion,  therefore,  of  a  small  amount  of  first-derivative  lead 
results  in  a  change  in  in  the  direction  of  more  lead.  A 

similar  analysis  of  the  effect  of  r  on  Sp  shows  that^;^=0 

when  r=0.  For  small  values  of  r,  therefore,  the  effect  of 
adding  first-derivative  lead  is  the  same  as  for  the  linear  servo: 
an  improvement  in  dynamic  stability  and  a  reduction  in  the 
amplitude  of  hunting  is  to  bo  expected.  For  large  values  of  r, 
however,  the  lag  due  to  increase  in  amplitude  (first  term) 
tends  to  counteract  the  direct  lead  (second  term)  and  may 
cause  more  violent  hunting.  Some  optimum  value  of  r  may 
lliereforo  exist.  A  direct  determination  of  its  value  may  be 
made  from  the  autopilot  for  various  vaku's  of  r,  6,  and  w. 

APPLICATfON  OF  METHOD  TO  CONTROLS  GEARED 
TO  MOTION  ABOUT  TWO  AXES 

Sometimes  tlie  aileron  or  rudder  is  geared  to  both  the  angle 
of  bank  and  the  angle  of  yaw.  Usually  this  gearing  is  ac¬ 
complished  by  tilting  the  gyro  in  such  a  way  that  the  rotation 
of  the  gyro  is  affected  by  bank  and  yaw  according  to  the 
formula 

4>  =  <f>  cos  T-{-ip  sin  r 


CONCLUDING  REMARKS 

The  frequency-response  method  of  analysis  is  a  useful 
graphical  means  of  determining  oseillation  eharacteristics  of 
an  airplane  equipped  with  an  automatic  pilot.  If  the  servo 
is  linear  (ideal  ease),  the  critieal  control  gearing  beyond 
which  increasing  oscillations  take  plaee  can  be  readily  de¬ 
termined.  If  the  servo  has  nonlinear  characteristics  (prac¬ 
tical  case),  the  existence,  amplitude,  and  frequency  of  steady 
hunting  oscillations  can  be  determined  approximately  from 
tbc  measured  frequency  response  of  the  servo  and  the  com¬ 
puted  frequency  response  of  the  airplane. 

In  general,  the  use  of  first-derivative  lead  (a  phase-shifting 
device)  has  a  favomble  effect  on  the  dynamic  stability, 
Large  amounts  of  first-derivative  lead  or  pha.se  shift  may  be 
destabilizing  due  to  the  increase  in  the  signal  amplitude  pro¬ 
duced  by  this  phase-shifting  device. 


Langley  Memorial  Aeronautical  Laboratory, 
National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  December  9,  1946. 


APPENDIX  A 


CALCULATION  OF  RESPONSE  OF  AN  AIRPLANE  TO  SINUSOIDAL  CONTROL  MOTION 


The  response  of  the  airplane  in  pitch  to  sinusoidal  elevator 
motion  will  be  used  as  an  example.  The  basic  equations  are 
derived  in  reference  3  but  with  somewhat  different  notation 
and  without  the  effects  of  control  deflection,  rate  of  change, 
of  angle  of  attack,  and  power.  The  equations  of  motion, 
for  a  given  control  deflection,  are 

{Xu-\-TD)u-\-Xaa-\-Xed=Q 

2-uU-\-  {zi — TD)d= — ziS  (Al) 

where 

_J_dT 

pVSdV 


^~pVS 


xe= 


O, 


2,1  —  Cl, 

{C.^+C,) 
_ Cl,  tan  To 

(y 

Cr.. 


c,„.~v 


dr,„ 


r 

'  in 


i>C,n 

da 


and 

To  trim  angle  of  climb,  deg 

T  thrust,  lb 

ly  pitching  moment  of  inertia,  slug-ft^ 

8 


m  mass  of  airplane,  slugs 

c  wing  chord,  ft 

S  wing  area,  sq  ft 

m=AF/F 

V  velocity,  ft/sec 

AF  change  in  velocity,  ft /see 
Cl  lift  coefficient 

a  angle  of  attack,  radians 

Cm  pitching-moment  coefficient 

C,y  drag  coefficient 

p  air  density,  slugs/cu  ft 

0 

Solving  equations  (Al)  for  -  by  the  method  of  determi¬ 
nants  gives 


tD 

0 

Za-\-xD 

—  2s 

d 

C,n^-\-Cmj^jD 

—  C,n^ 

1 

Xu-\-  tD 

Xa 

Xe 

2u 

Za-\-xD 

Ze—  tD 

C,n^ 

C,.^+C,n,D 

a.i,P-^yky^D^ 

Expanding  these  determinants  by  the  usual  methods  gives 
0 

an  expression  for  ^  consisting  of  the  ratio  of  two  polyno¬ 
mials  in  D,  which  may  be  written  as 

B (I  ]D~ h,T) C\ 

5  a-2D^2^boD^-\-C2D--\^d2D-\-e2 


The  phase  and  amplitude  of  (fie  response  of  6  to  the  motion 
5=sin  bit  is  obtained  by  substiluting  ibi  for  D  in  the  above 
0 

expression  for  This  substitution  gives  a  comple.x  number, 

say,  The  angle  of  lead  of  B  ahead  of  6  is  tan“* 

and  the  amplitude  of  6  is  that  is 

e,=  — tan  ‘  , 


and 


B=R  sin  {bit—e,) 


to 


APPENDIX  B 


DETERMINATION  OF  STABLE  SIDE  OF  STABILITY  BOUNDARY  OBTAINED  BY  FREQUENCY-RESPONSE  METHOD 


It  is  shown  in  the  body  of  this  report  that  a  critical  value 
of  the  control  gearing  exists  if  at  some  frequency  w  the  values 
of  €r  and  €p  are  equal  and  that  the  value  of  the  critical  control 
gearing  is  equal  to  at  that  frequency.  It  will  now  be 

shown  that  the  stable  side  of  the  stability  boundary  may  be 
determined  by  the  relative  slopes  of  the  €r  and  curves  at 
the  point  where  e,=ip. 

As  shown  in  reference  2,  the  response  of  an  airplane  to  any 
disturbance  may  be  computed  from  the  “response  function.” 
The  response  to  a  sinusoidal  disturbance  is  obtained  by  sub¬ 
stituting  iw  for  the  variable  D  in  the  response  function. 
The  response  functions  for  the  airplane  and  autopilot  may 
be  combined  to  form  the  stability  equation  of  the  airplane 
plus  autopilot  as  follows; 

Let  d=^fi(D)S  be  the  airplane  response  function  and 
S  =  kf2(D)6  be  the  autopilot  response  function  where  k  is  the 
control  gearing.  Combining  the  two  gives 

8  =  kMD)MD}8 

or 

\  —  kMD)fi{D)  =  Q 

By  snbstit.uting/(Z))=/i(I^)y2(f?)  this  equation  may  be  writ¬ 
ten  as 

l-kj{D)=Q  (Bl) 

This  is  the  “stability  ecpiation,”  the  roots  of  whieli  for  the 
variable  D  determine  the  dam])ing  and  frequency  of  the 
motion.  The  frequency  responses  are  obtained  from  the 
response  functions  by  use  of  the  following: 

fi  (iw)  =Re~''^=  \ 

Or 

f2{ioi)=8p^e^‘v 

Then  a  dynamic-stability  boundary  exists  if  1— ^/(ico)  =  0 
for  some,  value  of  w,  say  Wer-  It  follows  from  the  definition 
of  and/2  that  this  dynamic-stability  boundary  occurs  when 

er — 6^=0 

and 


at  some  value  of  u>=(x>cr-  Under  these  conditions,  equation 
(Bl)  has  a  pair  of  roots  D=  dzicocr  when  k  has  the  value  given. 
In  order  to  find  the  stable  side  of  the  boundary,  it  is  neces¬ 


sary  to  find  the  sign  of  the  I'eal  part  of  -^r  at  the  boundary. 


From  equation  (Bl) 


(ID  _  J{D) 
dk  kf{D) 


This  equation  must  bo  evaluated  when  D=ioicr-  At  this 
value  of  0),  f{io3cr)  =  }--  In  general, 

f^cr 


and 


If 


(litp—fr)' 
do) 


dD  . 


J  is  negative,  the  real  part  of  is  positive. 

As  k  increases  above  tlio  critical  value,  therefore,  the  roots 
of  the  stability  equation  indicate  instability;  that  is,  if 


<0 


then  the  stable  region  is  located  where  k  is  less  than  k^^  and 
vice  versa. 
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Positive  directions  of  axes  and  angles  (forces  and  moments)  are  shown  by  arrows 


Axis 

Moment  about  axis 

Angle 

Velocities 

Force 
(parallel 
to  axis) 
symbol 

Designation 

Sym¬ 

bol 

Designation 

Sym¬ 

bol 

Positive 

direction 

Designa¬ 

tion 

Sym¬ 

bol 

Linear 
(compo¬ 
nent  along 
axis)  : 

Angular 

Longitudinal 

X 

X 

Rolling  - 

L 

Y — 

Roll 

<i> 

0 

u 

P 

LatCTal  . 

Y 

Y 

M 

Z - *X 

V 
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Normal 

Z 

' 

Z 

N 

X - *Y 

w 

r 

Absolute  coefficients  of  moment 

/I _  ^  /I  _ ^  /I _  ^ 

'^1  «I.O  — _.cr  '^n — _IC 


ybS 

(roUmg) 


gcS 
(pitching) 


'qbS 
(yawing) 


Angle  of  set  of  control  surface  (relative  to  neutral 
position),  5.  (Indicate  surface  by  proper  subscrint.i 


4.  PROPELLER  SYMBOLS 


D 

P 

Diameter 

Geometric  pitch 

P 

P 

Power,  absolute  coefficient 

pID 

V' 

Pitch  ratio 

Inflow  velocity 

c. 

5  /  prs 

Speed-power  coefficient =.Wp^ 

V. 

Slipstream  velocity 

v 

Efficiency 

T 

Thrust,  absolute  coefficient  Ct= — Ssi 

n 

Kevolutions  per  second,  rps 

/  F  \ 

Q 

Torque,  absolute  coefficient 

Effective  helix  angle =tan~‘l2:j|^l 

6.  NUMERICAL  RELATIONS 


1  hp=76.04  kg-m/s=550  ft-lb/sec 
1  metric  horsepower= 0.9863  hp 
1  mph=0.4470  mps 
I  mp3=2.2369  mph 


1  lb=0.4536  kg 
1  kg=2.2046  lb 
1  mi= 1,609.35  m=5,280  ft 
1  m=3.2808  ft 
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ABSTRACT: 
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